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ABSTRACT 

The evolution leading to the formation of a neutron star in the very young Westerlund 1 star cluster 
is investigated. The turnoff mass has been estimated to be ^ 35 Mq, indicating a cluster age ~ 3 — 5 
Myr. The brightest X-ray source in the cluster, CXO J164710. 2-455216, is a slowly spinning (10 s) single 
neutron star and potentially a magnetar. Since this source was argued to be a member of the cluster, 
the neutron star progenitor must have been very massive (Mzams > 40 Mq) as noted by Muno et al. 
(2006). Since such massive stars are generally believed to form black holes (rather than neutron stars), 
the existence of this object poses a challenge for understanding massive star evolution. We point out 
while single star progenitors below M^ams ^ 20 Mq form neutron stars, binary evolution completely 
changes the progenitor mass range. In particular, we demonstrate that mass loss in Roche lobe overflow 
enables stars as massive as 50 — 80 M©, under favorable conditions, to form neutron stars. If the very 
high observed binary fraction of massive stars in Westerlund 1 (;> 70%) is considered, it is natural 
that CXO J164710. 2-455216 was formed in a binary which was disrupted in a supernova explosion such 
that it is now found as a single neutron star. Hence, the existence of a neutron star in a given stellar 
population does not necessarily place stringent constraints on progenitor mass when binary interactions 
are considered. It is concluded that the existence of a neutron star in Westerlund 1 cluster is fully 
consistent with the generally accepted framework of stellar evolution. 

Subject headings: stars: evolution - binaries: close - stars: neutron 



1. INTRODUCTION 

The recent detection of the variable X-ray source, CXO 
J164710. 2-455216, located 1.'6 from the core of the young 
open cluster, Westerlund 1, by Muno et al (2006) can 
be considered key to understanding the formation his- 
tory and nature of compact objects in young massive stel- 
lar clusters. CXO J164710.2-455216 is the brightest X- 
ray source in the cluster, radiating at a luminosity of 
~ 3 X 10^^ (D/bkpc)'^ ergs s~^ in the 0.5-8 keV energy 
band and characterized by a pulse period of 10.6 s. Based 
on subsequent XMM-Newton observations, Israel et al. 
(2007) and Muno et al. (2007) showed that the source 
increased its luminosity to ^ 10'^^ {D/5kpc)^ ergs s~^ and 
exhibited a rapid spin down after an outburst, exhibiting 
properties that are very reminiscent of highly magnetized 
neutron stars known as magnetars. 

Previous spectroscopic and photometric observations of 
Westerlund 1 were carried out by Clark et al. (2005), 
who showed that the population of massive stars could 
be understood if the cluster age was 4 ± 1 Myr (see also 
Crowther et al. 2006), a result which is also consistent with 
a more recent age determination based on an analysis of 
its intermediate and low mass stellar content (Brandner 
et al. 2008). For an assumed solar metallicity, a lower 
limit for the progenitor mass of CXO J164710. 2-455216 
was estimated to be ^ 35 M© . The reconstruction of the 
evolutionary history of CXO J164710. 2-455216 is highly 
desirable in placing limits on its progenitor, however, it is 
dependent upon whether it evolved in isolation as a single 



star or as a member of an interacting binary system. We 
note that the recent studies by Clark et al. (2008) indi- 
cate a high binary fraction, ^ 70%, for the massive Wolf 
Rayet star population in Westerlund 1 and that a binary 
evolutionary channel may be very likely. 

It is generally inferred from core collapse simulations 
that single stars more massive than ^20 — 25 Mq form 
black holes (e.g., Fryer 1999). Stars more massive than 
~ 20 — 25 Mq in binaries are also expected to form black 
holes in order to explain the Galactic black hole tran- 
sient systems and to satisfy nucleosynthetic constraints 
(e.g., Maeder 1992; Portegies Zwart, Verbunt & Ergma 
1997; Kobulnicky & Skillman 1997; Ergma & van den 
Heuvel 1998). However, a number of studies have been 
pointed out that binary interactions can significantly in- 
crease the neutron star/black hole mass formation limit. 
For example, Van den Heuvel & Habets (1984) studied the 
binary LMC X-3 containing a neutron star with a mas- 
sive companion, and found that stars as massive as 40 Mq 
may form neutron stars. Kaper et al. (1995) further in- 
creased this limit, estimating a mass of 50 Mq based on 
the analysis of a binary Wray 977, although it was later 
demonstrated that such a system could be explained with 
the much lower mass limit of 26 M© (Wellstein & Langer 
1999). However, Wellstein & Langer (1999) reemphasized 
that, in general, the critical initial mass limits for neu- 
tron star/black hole formation in binary systems may be 
quite different than for single stars and can be as high 
as 100 Mq. We point out that, in this regard, Brown et 



al. (2001) have noted that the final mass ol the compact 
object may depend upon whether the evolved core was 
exposed before or after central helium burning. More re- 
cently, Fryer et al. (2002) demonstrated that a star of 
60 M0 that was stripped of its envelope in a binary sys- 
tem can end up as a neutron star or a massive black hole, 
depending on the adopted wind mass loss rates during the 
Wolf-Rayet stage. These studies, taken as a whole, indi- 
cate that binary interactions can lead to a range of initial 
stellar masses separating the formation of neutron stars 
from black holes rather than a single limiting mass. We 
use this concept to examine the existence of a neutron star 
in the very young environment of Westerlund 1 . 

In this paper, we report on the formation history of neu- 
tron stars from very massive progenitors (^ 35 Mq) both 
as a single star and as a component of an interacting bi- 
nary system to determine the conditions under which a 
neutron star in a young cluster such as Westerlund 1 can 
form. In the next section, we outline the basic assumptions 
and input parameters of our numerical model as based on 
the StarTrack binary population synthesis code (Belczyn- 
ski et al. 2008). The numerical results are presented for 
the binary star and single star evolutionary channels in 
§3. The implications of these results for stellar and binary 
evolution of massive stars, as applied to Westerlund 1, are 
discussed in the last section. 

2. POPULATION SYNTHESIS MODEL 

Our population synthesis code, StarTrack, was ini- 
tially developed to study double compact object mergers in 
the context of 7 ray burst progenitors (Belczynski, Bulik 
& Rudak 2002b) and gravitational-wave inspiral sources 
(Belczynski, Kalogera, & Bulik 2002a). In recent years 
StarTrack has undergone major updates and revisions in 
the physical treatment of various binary evolution phases, 
and especially the mass transfer phases. The new ver- 
sion has already been tested and calibrated against ob- 
servations and detailed binary mass transfer calculations 
(Belczynski et al. 2008). It has been used in various appli- 
cations (e.g., Belczynski & Taam 2004; Belczynski et al. 
2004; Belczynski, Bufik & Ruiter 2005; Belczynski et al. 
2006; Belczynski et al. 2007). The physics updates that 
are most important for compact object formation and evo- 
lution include: a full numerical approach for the orbital 
evolution due to tidal interactions, calibrated using high 
mass X-ray binaries and open cluster observations, a de- 
tailed treatment of mass transfer episodes fully calibrated 
against detailed calculations with a stellar evolution code, 
updated stellar winds for massive stars, and the latest de- 
termination of the natal kick velocity distribution for neu- 
tron stars (Hobbs et al. 2005). For helium star evolution, 
which is of a crucial importance for the formation of dou- 
ble neutron star binaries (e.g., Ivanova et al. 2003; Dewi & 
Pols 2003), we have applied a treatment matching closely 
the results of detailed evolutionary calculations. If the he- 
lium star fills its Roche lobe, the systems are examined 
for the potential development of a dynamical instability, 
in which case they are evolved through a common enve- 
lope (CE) phase, otherwise a highly non-conservative mass 
transfer ensues. We treat CE events using the energy for- 
malism (Webbink 1984), where the binding energy of the 
envelope is determined from a set of He star models cal- 



culated with the detailed evolutionary code by Ivanova et 
al. (2003). In case the CE is initiated by a star crossing 
the Hertzsprung gap (HG) we assume a merger and abort 
further binary evolution. This is due to the fact that there 
is no clear core-envelope boundary (and no entropy jump 
as for more evolved stars) in the interior structure of HG 
donors to facilitate the formation of a remnant binary sys- 
tem. As a consequence, a large decrease in the formation 
efficiency of close double compact binaries results (Bel- 
czynski et al. 2007). Rejuvenation is considered in detail 
following the method presented by Tout et al. (1997). 
For a detailed description of the revised code we refer the 
reader to Belczynski et al. (2008). 

3. RESULTS 
3.1. Binary star progenitor 

We have calculated the evolution of 2 x 10^ massive bi- 
naries with our standard evolutionary model (Belczynski 
et al. 2008) for solar metallicity {Z = 0.02) and an initial 
mass function (IMF) with an exponent of —2.7 for pri- 
maries in the mass range 6 < Mzams,i < 150 Mq. The 
secondaries (4 < Mzams,2 < 150 M0) were chosen from a 
flat mass ratio {q = Mzams,2/-^zams,i) distribution. The 
initial separations were chosen from the distribution flat 
in logarithm (i.e., oc 1/a) limited by the maximum sepa- 
ration of 10^ Rq, while eccentricities were chosen from a 
thermal-equilibrium distribution (2e). 

As a result we have obtained populations of compact 
objects, i.e., white dwarfs, neutron stars and black holes, 
both single and in binaries. Single compact objects were 
formed upon disruption of a given binary during a super- 
nova explosion either through a natal kick, mass loss or 
both. We have used the kick distribution presented by 
Hobbs et al. (2005) represented by a single Maxwellian 
with a — 265 km s~^. We employed the above distribu- 
tion in its original form for neutron stars, while the kick 
velocity was decreased for black holes due to the fall back 
(e.g.. Fryer & Kalogera 2001). 

The stellar models employed in our population synthesis 
(Hurley, Pols & Tout 2000) result in an age of Westerlund 
1 of 5 Myr: this is a main sequence lifetime of a 35 M© 
star. Actual fits of single stellar models to over 50 massive 
stars in Westerlund 1 resulted in an age estimate of 4 ± 1 
Myr (e.g., Clark et al. 2005; Brandner et al. 2007). We 
adopt an age of 5 Myr, although we note that any estimate 
based on single star models is highly uncertain in the case 
of a cluster like Westerlund 1 with the very high content 
of binaries (^ 70%; Clark et al. 2008). At 5 Myr our cal- 
culations reveal a population of single neutron stars that 
originate from disrupted binaries. 

Evolution. In the following we briefly describe a typ- 
ical evolution leading to the formation of a very young 
solitary pulsar that originates from a disrupted binary. 
The evolution starts on the ZAMS with two massive stars: 
-^'-^zams.i = 65 Mq and Mzams,2 = 44 Mq, On a rather close 
a — 220 Rq and eccentric orbit: e = 0.65. The primary 
initiates Roche lobe overflow (RLOF) at t = 3.3 Myr while 
still on the main sequence (main sequence lifetime of an 
unaffected single star of that mass is Tms = 3.9 Myr, and 
it forms a He core of Mho = 18 Mq). At the onset of mass 
transfer (Case A) the orbit is circularized (a — 80 Rq , e = 
0) and the component masses (depleted by stellar winds) 



are Mi = 52 Mq and M2 = 41 Mq. The mass transfer 
is stable and proceeds on the nuclear timescale of the pri- 
mary with rates ranging from 2—8 xlO~^ M0 yr"-'^. We as- 
sume non-conservative evolution in which half of the mass 
lost by primary is accreted onto secondary, while the other 
half leaves the system, carrying away orbital angular mo- 
mentum. At the end of the primary main sequence phase 
we note the mass ratio has reversed (Mi — 31 M0 and 
M2 =^ 49 M0) and the orbit has expanded (a = 100 R©). 
The significant mass loss (RLOF) on the main sequence 
leads to the decrease of the central temperature, reducing 
the rate of hydrogen (H) burning in the core of the pri- 
mary. As a result, the primary forms a smaller core at the 
end of the main sequence phase (new Mhc = 10 Mq) and, 
thus, will form a lower mass compact object. Although the 
lifetime of the primary star increases with the decrease of 
its core mass, the fact that the mass transfer started late 
in the main sequence phase does not significantly affect 
the lifetime of the star (new Ty^s = 4.0 Myr). 

The mass transfer continues while the primary is cross- 
ing the HG (Case B), transferring mass at a high rate 
~ 10^^ Mq yr~^ such that the entire H-rich envelope is 
removed and the primary becomes a massive naked he- 
lium star with Mi = 10 Mq. The secondary, on the other 
hand, remains on the main sequence, but with its mass in- 
creased to M2 = 60 Mq. Mass transfer ceases as massive 
helium stars are characterized by small radii (i?i = 0.9 Rq 
for the primary). Thereafter, the primary evolves through 
the helium star phase, losing more mass via stellar winds 
and finally exploding in a Type lb supernova. At the 
time of explosion (t = A.SMyr) the primary has a mass 
of Ml = 6.6 Mq with core mass of 4.9 Mq. Such a star 
forms a massive neutron star (Mns = 1.9 Mq) as little 
matter is expected to fall back in such a case (e.g.. Fryer 
& Kalogera 2001). The explosion disrupts the system pro- 
vided that either a high natal kick or a kick off the orbital 
plane is applied, releasing the neutron star from the binary. 
The massive secondary soon (t = 6.1 Myr) completes its 
evolution and forms a solitary black hole of Mbh = 7.9 Mq . 

Physical properties. In Figure 1 we show the ini- 
tial mass distribution of very young pulsar progenitors. 
We note that the majority of the progenitors are massive 
-Wzams = 50 — 80 Mq sincc high initial masses are required 
in the framework of the above formation channel. Stars 
must be sufficiently massive to promptly complete their 
evolution within 5 Myr, but not too massive for otherwise 
black hole formation occurs. 

The mass distribution of the very young solitary pul- 
sars formed at iform <, 5 Myr is also illustrated in Fig- 
ure 1. As expected these pulsars are very massive: M^s = 
1.8 — 2.5 Mq. The high mass limit reflects our assumption 
on the maximum neutron star mass (Mns, max = 2.5 Mq), 
whereas the peak at Mns = 1-8 Mq originates from a wide 
range of progenitor masses. Specifically, our model for sin- 
gle stars has approximately a bimodal distribution of neu- 
tron stars masses; single stars with M^ams = 8 — 18 Mq 
form Mns = 1-3 Mq, while stars with M^ams = 18-21 Mq 
form Mns = 1.8 Mq. Such a distribution was adopted fol- 
lowing the study of Timmes, Woosley & Weaver (1996), 
which showed bimodality in the core mass of massive stars 
reflecting the differences in the entropy of their core re- 



gions. Therefore, since our scenario preferentially selects 
the most massive stars that form neutron stars the peak 
exists at Mns = 1.8 Mq, while there are virtually no neu- 
tron stars with Mns = 1-3 Mq. 

In Figure 2 the cumulative distribution of neutron star 
formation times is displayed with solitary neutron stars 
originating from binaries denoted by the solid line. Here, 
the fraction of neutron stars that have formed before a 
given time is shown. Neutron star formation begins at 
^form ^ 4.3 Myr, and is complete by tform ^ 60 Myr. The 
early phase of neutron star formation from very massive 
stars proceeded through the mass loss associated with the 
Roche lobe overflow episodes (as described above), while 
the neutron stars with the longest formation times re- 
sult from the evolution of intermediate mass stars either 
through electron capture supernovae (e.g. Podsiadlowski 
et al. 2004) or via rejuvenation. For the adopted age 
of Westerlund 1 i = 5 Myr we find that the fraction of 
neutron stars formed is /form = 0.001, implying that for 
1000 neutrons formed via disrupted binaries, one would 
form within the first 5 Myr. Since (i) most of the mas- 
sive stars in Westerlund 1 appear to reside in binaries 
and (a) most of the massive binaries (~ 95%; Belczynski, 
Lorimer & Ridlay 2008) are disrupted in supernovae explo- 
sions, one requires at least 1000 massive stars (^ 8 Mq) 
in the Westerlund 1 cluster to form one very young pul- 
sar. One should note the steepness of the cumulative curve 
for the times close to the age of Westerlund 1, where the 
fraction changes by more than two orders of magnitude 
from tform = 4.6 Myr (/form = 0.0001) to fform = 5.9 

(/form = 0.01). 

There are about ^ 10^ stars in Westerlund 1 and only 
stars with masses exceeding 2 Mq have had sufficient time 
to contract onto the main sequence. If we adopt a power- 
law IMF with an exponent of —2.7, there are about 10* 
stars more massive than ^ 8 Mq that can potentially form 
a neutron star. If these 10* massive stars are required to 
form the very young pulsar as observed in Westerlund 1 it 
would yield a fractional efficiency of 0.0001, corresponding 
to /form for a cluster age of 4.6 Myr. It should be noted 
that this approximate estimate is not intended to provide 
a calibrated rate/efficiency of young pulsar formation, but 
only to illustrate the an order of magnitude consistency of 
our results with the particular observation of one object. 

3.2. Single star progenitor 

The possibility that a primordial single star forms a very 
young pulsar has also been explored. The calculations 
were performed with the Star Track code for solar metal- 
licity and with standard wind mass loss rates as described 
in Belczynski et al. (2008). Within the framework of our 
assumptions regarding the relation between initial and fi- 
nal compact object masses (Belczynski et al. 2008), stars 
within the mass range of M^ams ^8 — 21 Mq can form 
neutron stars if the maximum neutron star mass is taken 
to be 2.5 Mq. The minimum formation time of a neutron 
star is tform = 9.4 Myr for the most massive progenitor 
(~ 21 Mq) and is longer by nearly a factor of 2 than the 
time required to produce a pulsar in Westerlund 1. This 
is shown in Figure 2 where the cumulative distribution 
of the neutron star formation times for primordial single 



'^AU details of the adopted scheme for neutron star mass calculation are given in Belczynski et al. (2008). 



stars (dashed line) is also shown. 

Adopting a higher maximum neutron star mass of 
3.0 Mq results in a slight decrease of the minimum for- 
mation time to tform = 9.0 Myr, corresponding to the fi- 
nal evolution of a Mzams — 21.5 M0 star. This is due to 
the fact that our predicted initial-final mass relation (see 
Fig. 1 of Belczynski et al. 2008) steeply rises with an initial 
mass above Mzams — 20 Mq as expected for the increased 
amount of fall back during compact object formation. 

A further decrease in the minimum time for formation 
of a neutron star could be accomplished by increasing 
the metallicity of the stars in Westerlund 1. For exam- 
ple, given Mns,max — 3.0 Mq and allowing for a metal- 
hcity oi Z = 1.5 x Zq = 0.03 as may be expected for a 
young cluster, the minimum time would decrease further 
to iform = 7.7 Myr (M^ams = 24.0 Mq), but still signifi- 
cantly high for the age of Westerlund 1. 

The only remaining alternative (within the framework 
of our adopted stellar models) is to increase the stellar 
wind mass loss rates. In this case, the mass loss rates may 
reduce the stellar mass sufficiently that even initially very 
massive stars form neutron stars rather than black holes. 
To examine this possibility, we have introduced two scal- 
ing factors to multiply our standard wind mass loss rates. 
Specifically, the factors /wind,H and /wind, Ho are used to 
scale winds of stars with H-rich envelopes and naked he- 
lium stars, respectively. In all the calculations presented 
herein, /wind,H = ./wind, He = 1- We follow the Hurley et 
al. (2000) scheme and the choice of wind mass loss rates. 
In particular, the rates from Kudritzki & Reimers (1978) 
and Iben & Renzini (1993) are used for giant branch and 
further evolution; Vassihadis & Wood (1993) are apphed 
to asymptotic giant branch; Nieuwenhuijzen & de Jager 
(1990) are used for massive stars, but modified for a metal- 
licity dependence (ex {Z/ZqY^'^); Hamman & Koesterke 
(1998) is employed for Wolf-Rayet-like winds; Humphreys 
& Davidson (1994) rates are used for luminous blue vari- 
ables (L > 6 X 10^ and R x \/L > 10'^; where L, R are 
luminosity and radius expressed in solar units). 

In Figure 3 the initial-final mass relation for single star 
models with various wind mass loss rates for stars with H- 
rich envelopes: /wind.H = 1, 2, 3 are presented. We have 
chosen models with high metallicity {Z — 0.03) and fixed 
wind mass loss rates for naked helium stars /wind. He = 1- 
Since the increased winds reduce the final compact object 
mass we note the change in mass of the most massive star 
that can still form neutron star. For Mns,nia.x = 3.0 Mq it 
is found that only stars with M^ams ^ 24, 40, 62 Mq 
can form neutron stars for /wind,H = 1, 2, 3, respec- 
tively. It is also found that only stars more massive than 
Mzams ^ 44, 50, 58 Mq can form compact objects within 
5 Myr. Both constraints can only be satisfied for the model 
with very high winds: /wind,H = 3. Therefore if single star 
is to produce a very young pulsar the stellar winds would 
need to be increased by factor of ~ 3 or higher. 

We find that such a drastic increase is unrealistic as 
it would prohibit formation of massive black holes over 
~ 4.5 Mq for a high metallicity and over ~ 6.5 Mq for 
solar metallicity. However, a few black holes in Galac- 
tic X-ray transient sources have dynamical mass estimates 
exceeding 10 Mq (e.g., Orosz 2003; Casares 2006). In ad- 
dition, recent observations indicate that stellar black holes 



can form with even higher masses in low metallicity envi- 
ronments; ^ 16 Mq black hole in M33 (Orosz et al. 2007) 
and ~ 23 M© black hole in ICIO (Prestwich et al. 2007; 
see also Silverman & Filipenko 2008 for an improved or- 
bital estimate). Although it is expected that black holes 
form with higher mass in lower metallicity environments, 
it was recently demonstrated that a reduction in the wind 
loss rate was necessary for the formation of such massive 
black holes (Bulik, Belczynski & Prestwich 2008). 

We have only investigated the increase of winds for stars 
with H-rich envelopes, as it has been argued that stel- 
lar wind mass loss rates for massive helium stars (i.e., 
Wolf-Rayet stars) are systematically overestimated due to 
the "dumpiness" of mass outfiows from these stars (e.g., 
Hamann & Koesterke 1998; Nugis & Lamers 2000). The 
reduction due to the dumpiness is already included in our 
study, but a further decrease (factor of ^ 2) was recently 
suggested by Vanbeveren, Van Bever & Belkus (2007), 
thereby, making it more unlikely that enhanced wind loss 
rates from single stars are responsible for the formation of 
neutron stars in Westerlund 1. 



4. DISCUSSION 

A population synthesis study of the massive star popu- 
lation for Westerlund 1 has been carried out to examine 
the possible formation history its neutron star member, 
CXO J164710.2-45526. The young age of the cluster (4±1 
Myr), corresponding to a star of 35 Mq at the turn off 
of the main sequence places important constraints on its 
prior history. It is found that an isolated neutron star can 
be formed as a result of the disruption of a binary system 
during a supernova phase (within 5 Myr) provided that 
the initial progenitor system is characterized by primary 
masses in the range of 50 — 80 Mq. On the other hand, 
only stars of mass ^ 60 Mq, for single star models can 
produce a neutron star within 5 Myr, however, the rate of 
mass loss via stellar winds must be enhanced to the extent 
that the formation of black holes with masses exceeding 
10 Mq becomes problematic. 

A major feature in the binary scenario of neutron star 
formation is the significant mass loss taking place in Case 
A RLOF. In this mass transfer phase, a much smaller He 
core (Mne ~ 10 Mq) is formed at the end of the main 
sequence phase as compared to the isolated evolution of 
the same star (Mne ^ 20 Mq). The removal of the H-rich 
envelope through the ongoing RLOF stage before central 
He-ignition (Case B) leads to the onset of strong Wolf- 
Rayet type winds, preventing an increase of the He core 
mass via H shell burning. Although we employ revised (de- 
creased to account for "dumpiness" ) wind mass loss rates 
for naked helium stars (e.g., Hamann & Koesterke 1998), 
the Wolf-Rayet phase leads to a further loss of ~ 3 Mq 
from the primary star. This is a secondary, but never- 
theless important effect leading to the formation of very 
young neutron stars. An additional factor enhancing the 
formation of neutron stars from very massive binary com- 
ponents was pointed out by Brown et al. (2001; and ref- 
erences within) , who argued that if the H-rich envelope is 
removed (e.g., in case B or early case C RLOF) and H shell 
burning is extinguished before the end of central He burn- 
ing the resulting Fe core mass is reduced in comparison to 
the case where the H-envdope is not removed. Hence, the 



formation of a low mass compact object (i.e., a neutron 
star rather than a black hole) could be enhanced. This 
finding was connected to the supply (or lack thereof in the 
case of the neutron star formation) of fresh helium into the 
burning core that can alter the amount of carbon and thus 
the carbon burning lifetimes in this phase of stellar evolu- 
tion. Since this effect is not modelled in our simulations, 
our predictions are more conservative. In other words, we 
could expect even more massive stars (than predicted here) 
to form neutron stars under the favorable conditions. 

Due to ejection of matter and natal kicks during a super- 
nova, it is expected that the neutron star will be displaced 
from its formation site. The projected distance of CXO 
J164710. 2-455216 from the cluster center is about 2.3 pc 
(e.g., Muno et al. 2006). The source could have been 
formed anywhere in Westerlund 1 and been moving since 
the time of its formation (iform) with velocity (Vns) im- 
posed on the neutron star during the supernova explosion. 
We note that motion of the system prior to the supernova 
explosion because of interactions with other stars or bina- 
ries in the cluster could also contribute, but the former are 
expected to be significantly larger. In Figure 4 the space 
velocities of the very young pulsars formed in our simula- 
tions under the assumption of a null initial velocity of the 
binary prior to explosion are shown. The results are pre- 
sented for two calculations; a reference calculation with the 
Hobbs et al. (2005) natal kick distribution {a — 265 km 
s~^), and one obtained with low kicks {a = 133 km s^^). 
It is found that the average velocities are Kis = 380, 210 
km s~^ for the reference and low kick models, respectively. 



In addition, due to the smaller probability of progenitor 
binary disruption (and thus formation of a solitary pul- 
sar) the number of young solitary pulsars decreases from 
-/Vns ~ 1200 for the reference model to ~ 800 for the low 
kick model. Using the velocity distribution with the age 
of each pulsar (Fig. 2), the average distance traveled since 
its birth is calculated as D^s = 90, 50 pc for the reference 
and low kick models, respectively. This reveals that the 
majority of the young pulsars is bound to leave the clus- 
ter. For a region characterized by 10 pc (comparable to 
the size of Westerlund 1) an estimate of the probability of 
finding a pulsar within the cluster yields ~ 10%, 20% of 
pulsars for the reference and low kick models, respectively. 
This result may indicate that the predicted fractional ef- 
ficiency of forming very young pulsars is lower (by factors 
of '-^ 5—10) than presented in Figure 2. The decreased 
efficiency (even by factor of 10), however, still permits one 
very young pulsar to be found within the cluster age of 5 
Myr. 
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Fig. 1. — Top panel: Initial (Zero Age Main Sequence) mass of neutron star progenitors as predicted in population 
synthesis calculation for Wester lund 1 at the age of 5 Myr (solar metallicity, standard winds). All neutron star progenitors 
originate from binaries that were disrupted in a supernova explosion. Note the very high masses of the progenitors 
{M^ums ~ 50 — 80 M0). Bottom panel: Neutron star masses predicted in the same calculation. Note the high mass 
{Mns ^ 1-8 — 2.5 Mq) of the neutron stars. 
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Fig. 2. — Cumulative distribution of formation time of single neutron stars predicted in population synthesis calculation 
for Westerlund 1. Single neutron stars that originate from disrupted binaries are shown with solid/blue line, while the 
ones that originate from primordial single stars are shown with dashed/red line. Note that primordial single stars do not 
form neutron stars below 9 Myr, while a small fraction of neutron stars (/form ~ 0.001) form in the first 5 Myr from 
binary progenitors. 
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Fig. 3. — Initial-final mass relation for single star models with various wind mass loss rates for stars with H-rich envelopes: 
/wind,H = 1, 2, 3. All models are presented for high metallicity {Z = 0.03) and fixed wind mass loss rates for naked 
helium stars (/wind, He = !)• We mark the initial mass (Mzams) range for each curve (shaded) for which formation of a 
compact object takes shorter than tform — 5 Myr. Also marked is the final compact object mass (Mf) range in which it 
is expected that the transition of a neutron star to black hole formation takes place. Note that only for models with the 
very high winds (/wind,H = 3; bottom curve) can neutron stars be formed from single progenitors within 5 Myr. 
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Fig. 4. — Neutron star space velocities for two natal kick models; solid line shows results of our reference calculation with 
Hobbs et al. (2005) kicks (cr = 265 km s^^), while dashed line represents results obtained with low kicks (cr = 133 km 
s"^). Note that not only the average velocity decreases from reference (Kis ~ 380 km s~^) to low kick model (1/ns ~ 210 
km s~^), but also the number of young solitary pulsars (A^ns '^ 1200) decreases due to the smaller probability of progenitor 
binary disruption for low kicks (A^ns ~ 800). 



